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Abstract. Thermodynamic framework of finite strain viscoelasticity with 
second order weak nonlocality in the deformation gradient is investigated. The 
application of Liu procedure leads to a class of third grade elastic materials 
where the second gradient of the stress appears in the elastic constitutive rela- 
tion. Finally the dispersion relation of longitudinal plane waves is calculated 
in isotropic materials. 



1. Introduction 

Thermodynamic requirements are important in all theoretical approaches of con- 
tinuum mechanics. The classical form of the Clausius-Duhem inequality [I] does 
not allow higher than first grade elasticity [2] and plasticity [3] without any further 
ado. Therefore, new concepts - like higher order stresses, configurational forces - 
emerged to circumvent this condition and understand the thermodynamic compat- 
ibility of successful material models [H |B] . 

In this paper we show that a weakly nonlocal extension of the constitutive state 
space does not contradict to the Second Law and leads to constitutive relations 
of higher grade finite strain elasticity and viscoelasticity that are compatible with 
rigorous thermodynamic methods and requirements. Our method is based on two 
basic observations, that are different from the classical framework of Gurtin [5] 

- the entropy flux is a constitutive quantity, 

- for higher order weakly nonlocal state spaces the gradient of the balances 
and other kinematic constraints result in further constraints on the entropy 
inequality, 

The assumption of a constitutive entropy flux is a straightforward generalization 
of the Gibbs-Duhem inequality and one can prove that in simple cases leads to 
the classical form and to the classical results both in irreversible thermodynamics 
[7] or in thermoelasticity in particular [5] . This generalization is well accepted and 
applied beyond mechanics [9j[T0j[Tl]. With the assumption of a constitutive entropy 
flux the Liu and Coleman-Noll procedures are equivalent (see [12] for a proof in a 
particular case). 

It is also remarkable that further constraints in the Liu procedure result in 
more general constitutive functions 13 . Our main result, the thermodynamic 
admissibility of the dependence of constitutive functions on space derivatives of the 
deformation gradient, is the consequence of this general property of the entropy 
inequality. Therefore we do not need to introduce higher order stresses (neither 
any other additional physical concepts) in advance, but at the end we will see, that 
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some of the consequences of our method can be interpreted in that terms. The 
introduction of the gradient of local constraints (e.g. balances) in non-equilibrium 
thermodynamics is a mathematical necessity in case of higher order weakly nonlocal 
constitutive state spaces - overlooked by Gurtin in [2] - and has a unifying power 
to understand the role of the Second Law in several seemingly different theories of 
continuum physics (see e.g. |14j and the references therein). 

In this paper we apply our method with Liu procedure in third grade elasticity 
(DP i p63), where the constitutive state space depends on the second space derivative 
of the deformation gradient. Therefore third grade elasticity is classified as a second 
order weakly nonlocal theory. We prove the thermodynamic admissibility of a class 
of constitutive relations with two remarkable properties 

- the second order derivatives of the deformation appear without explicitly 
introducing double stress in advance as an independent theoretical concept, 

- the second derivative of the stress is part of the nondissipative stress-strain 
constitutive relation. This is similar to the suggestion of Aifantis (see [T5] 
and the references therein). 

We also derive a simple dispersion relation of longitudinal plane waves to demon- 
strate the properties of the constitutive relation. 

2. Continuum in a Piola-Kirchhoff framework 

All quantities are defined on the reference configuration. The substantial time 
derivative is denoted by a dot and the material space derivative is di, where i 6 
{1, 2, 3}. Higher order derivatives are denoted by more indices, e.g. dij is the second 
gradient, x 1 is the motion, — dj\ l is the deformation gradient. 

According to the traditional concept of objectivity |16j this kinematic standpoint 
ensures the objectivity of the whole treatment as long as the constitutive quanti- 
ties depend on objective physical quantities. However, the original mathematical 
formulation of the concept of objectivity by Noll is questionable [HI [18] and a gen- 
eralization based on precise spacetime notions and a four-dimensional formalism 
was suggested to improve it [HJ [2UJ . In this work we do not apply this generalised 
objective framework, our results are derived by usual three dimensional notions. 
However, we exploit some consequences of this generalization to simplify our cal- 
culation. First of all it will be convenient to work in a Piola-Kirchhoff framework, 
where the balances and the physical quantities are interpreted in the reference con- 
figuration (see [U(2I] for a similar treatment). The first Piola-Kirchhoff stress will 
be denoted as a tensor. As we have mentioned above, a constitutive state space 
spanned by Noll-objective physical quantities (e.g. right Cauchy-Green deforma- 
tion) the objectivity of the whole treatment could be ensured. However, it is more 
convenient to work with the deformation gradient, the material velocity and the 
total energy as basic variables. Moreover, it is not forbidden by our generalized 
notion of objectivity. We will partially change to Noll-objective quantities at the 
end introducing the internal energy with the usual definition and show a particular 
stress-strain relation with a Cauchy deformation measure. 

Therefore in our treatment of the constitutive theory of third grade elastic 
materials the constitutive state space is based on the following fields (v l ,djV l , 
djkv 1 , F l j,dkF l j, dkiF l j,e, die). Here v % is the velocity, F\ is the deformation gradi- 
ent, and e is the specific total energy. Our approach is second order weakly nonlocal 
in the velocity and the deformation gradient, and first order weakly nonlocal in the 
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energy. It is usual to avoid introducing velocity field by working with internal en- 
ergy and internal energy balance as an independent variable. However, the velocity 
and deformation gradient fields form a single physical quantity therefore we find 
instructive to show that the direct definition of internal energy can be introduced 
at the end, and that the total energy filtered through Liu procedure can give the 
same results when specified to local constitutive relations. This is the approach that 
we have followed in case of relativistic fluids, where Liu procedure was essential to 
distinguish total and internal energies [35]. 

The well known kinematic relation between the velocity and the deformation 
gradient 

(1) F) - djV* = 

is introduced as a constraint for the entropy inequality together with the balances. 
The balance of linear momentum is 

(2) /<,;./•' - (1,1" = ()'. 

Here po is the material density and T 1 ^ is the first Piola-Kirchhoff stress, introduced 
as a tensor. The balance of total energy is 

(3) Po e + dtf = 0, 

where q l is the energy flux. As we are working in a second order weakly nonlocal 
constitutive state space the derivative of the kinematic relation ([T]) and the momen- 
tum balance ((2]) are additional constraints according to the exploitation method of 
weakly nonlocal continuum theories [23l [24] : 

(4) <w; ; o,y = o%, 

(5) podjv* + d jk T ik = I)',. 

The gradient of the energy balance does not give an additional constraint, be- 
cause the constitutive state space is first order weakly nonlocal in the energy. 
The entropy inequality requires that 

(6) Po s + diJ 1 > 0, 

where s is the specific entropy and J 1 is the material entropy flux. Here we are look- 
ing for restrictions on the constitutive functions T lJ ,q 2 , J 1 in terms of the entropy 
density derivatives. It is important to see, that the derivative of the momentum 
balance extends the process direction space, which is spanned by the first and also 
the second space derivatives of the constitutive state space. 

We introduce A|, Aj, K, A^ fe , \{ Lagrange- Farkas multipliers of the equations (p}- 
([5]) respectively. 
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Therefore the starting point of the Liu procedure is the following inequality 

< Po s + d t J l - A{{F) - djV*) - Xiipov* - fyTV) - K ( Po e + dtf)- 
Ai k (8 k F) - d kj v l ) - Xiipofyv* + d k] T lk ) = 



Po 



ds ds ds ds * ds 

' v l + po^r-jdji 1 + Pa — — -d jk v l + Po^—F) + Po ^ d k F\ + 



'dd k F\ 



<>(- ddjv uw JK . . , ; , 

ds „ ■. ds . ds 



'OF) ; 



Po 



dd k ,F 



-d k iF) + PoT-e + pr ) —-d i e+ 
J rip. rirhp 



d d J k BP d J k d 7' 

dJ m - „ dJ* dP 
ddi k F\ 3 de dd t e 



J ■ dq 3 i 



1 d kl v l + J^-d m v l + ^-ih /• ', + F\ - 



up 

dq m 
dd lk F 



ddkjV* 

dq^ 

3 de 



dd k F) 



-dmikF^ + —d t e + -^-djie 



/ dT rj 



dT 



rk 



ddjV' 



■d kj v l 



ddi k F' 



-dmlkF) 



dT rl 

dd kj v 
dT 



■dikjv 1 



dT 



rk 



dF) 



■d k F) 



dT 



rl 



dd k F) 



-d lk F)- 



de 



die - -——d^e 



ddie 



A{(F)-dyy 



K ( PodsV 1 

Q T rl 



dd kj v i 



dv 

dsikjv 1 - di kj v l d; 



dv 1 



dT 



r k 



ddjV 1 



d skj v l - d k ,v l d. 



kj 1 



dT 



rk 



ddjV 1 



dT rl 

dd ki v 



dF' 



-d sk F) - d k F)d s 



dT 



rk 



dF 



dT rl pi 
dd k F) slk 3 



d lk F)d s 



dT 



rl 



dd k F\ 



dd lk F 



-dsmlkF^ - dmlkF^dg 



3 

Qrprm 

dd lk F\ 



de 



■d„. 



die d„ 



dT r 



Q T rj 

-— - d s jie dj^e d s 
dd,e 



(7) 



de 

Ai k (d k F)-d kj v 



Q T rj 

ddie 



The time derivative related Liu equations imply that Lagrange-Farkas multipliers 
are determined by the entropy derivatives, and that the specific entropy and the 
stress do not depend on the highest derivatives of the constitutive state space 
s = s{v\d 3 v\F),d k F),e). 

Then, the Liu equations related to the highest order space derivatives determine 
the entropy flux in the following form 



( 8 ) Jt -P-i§h 

de dd t v r 



i>I " ., dT kr ., _ dT kr „ , ,, 



de 



ddiv' 
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Here K l — K z (v l ,djV l ,F l : j,dkF l j ) e) is the extra entropy flux with the above de- 
noted restricted functional dependencies. Then we introduce several convenient 
assumptions in order to get a simple and solvable form of the dissipation inequal- 
ity. First of all we define the internal energy of the third grade viscoleastic material 
including an isotropic kinetic energy contribution of deformation gradient rate and 
also the classical heat flux q: 



("I ^=e-^i-y(i^) 2 -f ^F{, <}<: -.ji + r.T". 



Morever, we may observe that a particular choice of the extra entropy flux reduces 
the dissipation inequality to a solvable form. Therefore we assume that 

K* := %^ - -^ T d j v k . 
de J ddiF k . 3 

Finally the temperature 6 is defined by the entropy derivatives §§ = In = g an d 
the free energy as ^>(F, VF) := u — 9s. Then we obtain 

(10) 

0a s = 6q%l + drf (t? - || - a ldlk T l H{ - a 2 d\T k t + dk^A > 0- 

The first two terms in the parenthesis of the above expression are the classical 
terms from second grade elasticity. The very last term in (fTU)) resembles the double 
stress relation that one can get by virtual power techniques. 

In the non dissipative case assuming constant material parameters we get a 
constitutive relation of third grade elasticity in the following form 

(If) T-V-(a 1 (V-T)I + a 2 VT) = ||-V. 



OF <9VF ' 



3. Simple waves 

We may calculate the dispersion relation of a one dimensional plane wave in 
the small strain approximation in order to check some consequences of the above 
stress-strain relation. 

Let us assume a usual isotropic quadratic free energy for symmetric strains e* = 

5 j + — ) ' *^ a * * s secon d order isotropic quadratic also in the gradient of the 
symmetric strains. Then according to representation theorems only two additional 
material parameters a\ , a 2 appear in the free energy function [25] 

(12) ^),d k e)) = ±(tf+^e) + ^di^&el + ^d^d'e). 
Then we get 

(13) Tl - aMT^Sf - a 2 d k &T\ = \f%6> + 2 M e| - a^DH ~ a 2 d k (4). 

Let us investigate the simplest one dimensional case and reduce the treatment to 
one component of the above tensorial equation. Introducing the notation T = T 11 
and e = e 11 and also denoting d\ by a dash we get 



— al = Ae — ae . 
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Figure 1. The square of the phase velocity v 2 = (u)(k)/k) 2 with 
the parameter values po — 1 and A = 1, a = 1 and a = 2. 

where a — ct\ + a%, X = X + 2p and a = a\ + &2- This constitutive relation is 
coupled to the balance of momentum that is in our case 

p e - T" = 0. 

Assuming constant material parameters we get the following dispersion relation 

2 _ k 2 (X + ak 2 ) 
po(l + ak 2 ) 

The characteristic feature of this dispersion relation is that the small wavelength 
and large wavelength limits result in finite acoustic phase velocities: 

linifc^o — J X/po and limfe_ i . 00 = y / a/(apo) as it is demonstrated on 
Figure [TJ This kind of behaviour is a property of the double wave equation [26) . 
Double wave equations are introduced by microstructural considerations e.g. in 
microstrain theories or internal variable theories [551 123 ■ 

4. Summary 

In this paper we have investigated a weakly nonlocal extension of viscoelasticity 
up to second order in the deformation gradient. The entropy flux was considered as 
a constitutive quantity and we applied Liu procedure introducing also the gradients 
of the balance of momentum and of the kinematic relation ([lj as constraints of the 
entropy balance. 

The calculations were performed in the Piola-Kirchhoff framework. The consti- 
tutive state space was chosen according to the generalization of the Noll principle 
of frame indifference introducing the components of the velocity-deformation gra- 
dient mixed four-tensor and their first and second space derivatives as constitutive 
variables. 

We have given a complete solution of the Liu equations and obtained a partic- 
ular form of the entropy flux and the dissipation inequality. In order to solve the 



THERMODYNAMIC CONSISTENCY OF THIRD GRADE FINITE STRAIN ELASTICITY 7 



dissipation inequality we have introduced some simplifications. In particular a qua- 
dratic form of the kinetic energy and some further related rearrangements resulted 
in a form where Onsagerian fluxes and forces can be identified. The constitutive 
relation of the non-dissipative mechanical material contains the gradients of the 
pressure in addition to the classical terms and the usual form of the constitutive 
relation of the double stress appeared without postulating such term in advance. 

It is important to note that the extension of the constitutive state space toward 
higher order gradients does not change the obtained form of (|10p as long as only 
the first derivatives of ([1]) and @ are introduced as additional constraints. This 
extension results in simple explicit solutions of the non-dissipative differential stress- 
deformation relation and then can be considered as a constitutive relation of 
third grade elasticity. 

We have calculated also a one dimensional dispersion relation and concluded, 
that it is similar to the dispersion relation of some microstructured materials as 
one could expect in case of higher grade solids (see e.g. |28|). 

As we have mentioned in the introduction, stress-strain relations similar to (| 1 1|) 
were already proposed in the literature |15j . However, there the motivation was to 
remove stress singularities with a kind of ad hoc "reaction diffusion" form. Here 
we have shown that this kind of extension is compatible with a weakly nonlocal 
thermodynamic framework, there are natural boundary conditions coming from the 
requirement of vanishing entropy flux and a physical origin of stress derivative terms 
is the modified (isotropic) kinetic energy, related to the rate of the deformation 
gradient. 
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